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ABSTRACT 

We perform a statistical weak lensing or "cosmic shear" analysis of the COMBO- 1 7 survey - a 
unique dataset with shear quality R-band imaging and accurate photometric redshift estimates 
(az — 0.05) for ^ 90% of galaxies to ttir < 24.0. We undertake a full maximum likelihood 
analysis to measure directly from the data the weak lensing power spectra, C^", and 

C^'^ in 5 band powers from £ = 400 to £ = lO'^, where K is the usual lens convergence 
and P is an odd-parity 'curl' component of the shear signal. We find a strong measurement 
of the convergence power over five fields. The non-gravitational /3-field has a much lower 
significance, indicating our data is free of major systematics, while the cross-correlation of 
ti and (3 is consistent with zero. We have also calculated the shear correlation functions and 
variance over a range of scales between 0.5 and 20 arcmin. Our measurements of minimal 
star-galaxy correlations and cross-correlations between galaxy components provide further 
evidence that any systematics are negligible. 

In addition, we have used our results to measure cosmological parameters, constraining 
the normalisation of the matter power spectrum to be as — (0.72 ± 0.09)(rim/0.3)~" '*^, 
where the errors quoted are 1-ct due to the intrinsic dispersion in galaxy ellipticities, cosmic 
and sampling variance. We have significantly reduced the usual additional uncertainty in the 
median redshift (z„j) of the source galaxies by estimating z„i directly from our data using ac- 
curate photometric redshift information from the COMBO-17 multi-band wide-field survey. 
To demonstrate the power of accurate redshift information, we have also measured parameters 
from a shear analysis of only those galaxies for which accurate redshift estimates are avail- 
able. In this case, we have eliminated the uncertainty in the redshift distribution of sources 
and we show that the uncertainty in the resulting parameter constraints are reduced by more 
than a factor of 2 compared to the typical uncertainties found in cosmic shear surveys to date. 
Finally, we combine our parameter measurements with constraints from the 2dF Galaxy Red- 
shift Survey and with those from the CMB. With these additional constraints, we measure the 
normalisation of the matter power spectrum to be as = 0.731qq3 and the matter density of 
the Universe to be 0„, = 0.27to:E!?- 

Key words: cosmology: observations - gravitational lensing, large-scale structure of the Uni- 
verse 



1 INTRODUCTION 

There continues to be great interest and significant progress in 
measuring the weak lensing signal arising from large-scale struc- 
ture. This phenomenon, observed as the weak coherent distortion 
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of background galaxies due to light ray deflection by interven- 
ing matter, offers us a direct probe of the mass distribution in the 
Universe. We can consequently measure cosmologically important 
quantities (see e.g.. Bemardeau et al. 1997; Jain & Seljak 1997; 
Kamionkowski et al. 1998; Kaiser 1998; Hu & Tegmark 1999) such 
as the bias (see e.g. Hoekstra et al. 2001) and the normalisation of 
the matter power spectrum (e.g. Bacon et al. 2002; Hoekstra et al. 
2002; Refregier, Rhodes & Groth 2002; Van Waerbeke et al. 2002). 
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Of particular interest is the prospect of reconstructing the 
power spectrum of mass fluctuations from the weak lensing signal. 
Pen et al. (2002) have used an estimator for the weak lensing shear 
power spectrum, obtainable from correlation function measure- 
ments, to measure the shear power spectrum from the DESCART 
survey, while Schneider et al. (2002) have further developed es- 
timators for this purpose. Meanwhile, Hu & White (2001) have 
demonstrated the utility of a maximum likelihood approach to re- 
constructing the shear power spectrum. Here we will apply a full 
maximum likelihood analysis to a cosmic shear survey for the first 
time. 

In this paper, we describe the results of a weak shear analysis 
of the COMBO- 17 dataset, acquired with the La Silla 2.2m tele- 
scope in Chile (Wolf et al. 2001). This survey includes 1.25 square 
degrees of deep R-band observations from which we draw our sam- 
ple of background galaxies for shear measurements. In addition, the 
survey has yielded photometric redshifts for approximately 40% of 
the objects in our selected sample of galaxies, which we use to im- 
prove our understanding of the shear signal. We apply a maximum 
likelihood reconstruction to our dataset to obtain the shear power 
spectrum, and we fit cosmological models to measurements of both 
the shear power spectrum and the correlation functions to measure 
joint constraints on the matter density Q.m and normalisation of 
the matter power spectrum ug- We obtain further constraints on 
these parameters by combining our measurements with constraints 
from both the 2dF Galaxy Redshift Survey (2dFGRS, Percival et 
al. 2002) and from CMB data (e.g. Lewis & Bridle, 2002). 

Cosmic shear studies require careful data reduction followed 
by assessment and removal of systematic effects, such as shear 
induced by a telescope, anisotropic point spread function (PSF), 
and circularization of galaxy shapes. We demonstrate in this paper 
that our shear catalogues have almost negligible systematic errors 
after applying appropriate corrections, and we carefully audit the 
remaining sources of error (from e.g. sample size, redshift uncer- 
tainty and shot noise) in order to interpret our results. 

The paper is organized in the following fashion. In Section 2 
we summarize the necessary formalism for our weak lensing anal- 
ysis, including definitions of shear and the shear power spectrum. 
In Section 3 we describe the COMBO- 17 survey. We then discuss 
the procedures used for data reduction, and describe our method 
for obtaining a shear catalogue fully corrected for systematic ef- 
fects. Section 4 examines the individual COMBO- 17 fields in de- 
tail, and explains our approach to including redshift information for 
these fields. We measure shear correlation functions and cell vari- 
ance in Section 5, checking for residual systematics and comparing 
with the findings of other groups. Section 6 contains our likelihood 
analysis, resulting in shear power spectra measured for single fields 
and for the survey as a whole. In Section 7 we use our correlation 
functions and power spectra measurements to estimate cosmologi- 
cal parameters, taking into account the errors from sample size and 
redshift uncertainties. We then combine our constraints on these pa- 
rameters with constraints obtained from the 2dFGRS and the CMB. 
In Section 8 we discuss these results and summarize our conclu- 
sions. 



2 WEAK LENSING QUANTITIES 

In the last few years weak gravitational lensing has emerged as the 
most direct method for measuring the distribution of matter, regard- 
less of its nature, in the Local Universe (Mellier 1999, Bartelmann 
& Schneider 2001). This is largely due to its basis in well under- 



stood physics: weak lensing is essentially a scattering experiment 
of photons off the gravitational field generated by cosmological 
structure. Here, we describe the various weak lensing fields and 
their statistical properties. 



2.1 Weak lensing fields 

Weak gravitational lensing induces a distortion into the images of 
distant source galaxies. This distortion can be parameterised by 
measuring the ellipticities. 
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of each galaxy from its quadrupole moments at a given isophotal 
threshold (Kaiser, Squires and Broadhurst, 1995) or by examining 
the distortion to a set of orthogonal modes describing the galaxy 
shape (Refregier & Bacon 2003; Bernstein & Jarvis, 2002). In the 
first instance - which is the approach we adopt in the following 
analysis - the effect of lensing is to induce an additional ellipticity 
on the galaxy image; 



= eij + 27i. 



where 
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(3) 



is the trace-free lensing shear matrix. The components 71 and 72 of 
the shear matrix represent the two orthogonal modes of the distor- 
tion. The shear matrix is recoverable since in the absence of intrin- 
sic alignments of galaxies, the galaxy ellipticities average to zero, 
(e) = 0. 

Since gravity is a potential theory in the weak-field regime, 
the shear field can be related to a lensing potential 



la = i did J - 



(4) 



where di = r{5ij — fifj)Vi is a dimensionless, transverse dif- 
ferential operator, and = did' is the transverse Laplacian. The 
indices (t, j) each take the values (1,2), and we have assumed a flat 
sky. On the scales currently of interest (i.e. from ~ 100 h~^'k\)c to 
100 /i~^Mpc), this is an excellent approximation. 

The lensing potential is also observable via the lens conver- 
gence field 



(5) 



which can be estimated from the weak magnification of sources. 



1 + 2k. 



(6) 



This can be measured from either the change in galaxy number 
density (Broadhurst, Taylor & Peacock 1995; Taylor et al. 1998) or 
from the change in size of sources (Bartehnann & Narayan, 1995). 
Although this has a lower signal-to-noise than the shear distribu- 
tion, it is an independent estimator of the lensing potential. 

The convergence field is related to the shear field by the dif- 
ferential relation, first used by Kaiser & Squires (1993); 



d didj'fij, 



(7) 



where d ^ is the inverse 2-D Laplacian operator defined by 

d~'^ j (fr ln\f - r\. (8) 
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A useful quantity for tracing noise and systematics in gravita- 
tional lensing is the divergence-free field, f3, defined by 



(3 = d '^Eidjdn-yij. 

where e" is the Levi-Civita symbol in two dimensions, 



-1 

1 



(9) 



(10) 



If 7ij is generated purely by the lensing potential, this quantity 
vanishes. But if there are non-gravitational sources, due to noise, 
systematics or intrinsic ahgnments, f3 will be non-zero. In addition 
/3 terms can arise from finite fields, due to mode-mixing of the k 
and 13 fields (Bunn, 2002). 

In a spatially flat Universe, with comoving distance r we can 
relate the lensing potential and the gravitational potential by (e.g. 
Kaiser 1998; Hu 2000) 



(r) = 2 / dr' 



rr' 



(11) 



This equation assumes the Bom approximation, in which the path 
of integration is unperturbed by the lens. 

Although the lensing potential depends on the distance to the 
source galaxy, this dependence is usually lost by averaging all lens- 
ing quantities over the source distribution. However with redshift 
information the full 3-D character of lensing can be usefully recov- 
ered (Taylor 2001, Bacon & Taylor 2002, Hu & Kccton 2002). 

Finally the Newtonian potential, <I>, can be related to perturba- 
tions in the matter density field, S = Spm/ Pm, by Poisson's equa- 
tion; 



(12) 



where we have introduced the cosmological scale factor o, the Hub- 
ble length \h = c/Ho « 3000 /i~^Mpc, and the present-day 
mass-density parameter Qrn- Here, Ho is the Hubble constant and 
c is the speed of light. 



2.2 Statistical properties 

2.2.1 Shear covariance matrix 

We may define a shear covariance matrix by 



(13) 



where the indices (a, h) each take the values (1,2). Fourier trans- 
forming the shear field. 



(14) 



and decomposing it using equations (7) and (9), we may generate 
the shear power spectra from correlations of n and (3: 

{K{t)K{l')) = {2TTfCr5D{£-£'), 

{n{l)p*(l')) = {2^fcf5D(l-l'). (15) 

The parity invariance of weak lensing suggests that C^'' = C^'' = 
0. However other effects, such as noise and systematics, as well 
as intrinsic galaxy alignments may give rise to a non-zero C^^ . 
Hence in our analysis we shall leave its amplitude to be determined 
by the data. The cross-correlation of «;(£) and /3(£) is expected to be 
zero but it will also allow a second check on noise and systematics 



in the shear field, and we shall treat is as another free function. 
In particular finite field and boundary effects can lead to leakage 
of power between these three spectra, which we shall attempt to 
monitor. 

The shear power spectrum and the convergence power are re- 
lated by 



(16) 



in the flat-sky approximation. For a spatially flat Universe, these 
are in turn related to the matter power spectrum, Ps{k, r) by the 
integral relation (see e.g. Bartelmann & Schneider 2001): 



W{r) 



(r) 



(17) 



where a is the expansion factor and r is comoving distance, rn is 
the comoving distance to the horizon: 



rn = c 



dz 

wry 



(18) 



where the Hubble parameter is given in terms of the matter density, 
fim, the vacuum energy density, Q.v and the spatial curvature, 

as 



H{z) = Hq[{1 + zfn„, + (1 + zf^K + 



il/2 



(19) 



The weighting, W, is given in terms of the normalised source dis- 
tribution, G{r)dr = p{z)dz: 



W{r)= I dr' G{r') ''—^. 



(20) 



The covariance of the components of the shear field are related 
to the power spectra by (Hu & White 2001) 

C7ii(f) = J-0^{crcos^'i<Pe + C^^sm^^^i 

+Cf 8m4ife)\W(£)\'e'^-^, 

+Cfcos4^e)\W{l)fe'^'^, (21) 

where cos ipe — £.£x and £,r is a fiducial wavenumber projected 
along the x-axis. We have included here a smoothing or pixelisation 
window function; 



W{i) = jo(46lpix/2)jo(4^pix/2), 



(22) 



where jo = sin(x) /x is the zeroth order spherical Bessel function 
and 0pix is the smoothing/pixel scale. 



2.2.2 The rotated shear correlation function 

Another important two-point statistical measure of the weak lens- 
ing signal are the correlation functions, defined by 



C2{e) = (72(f)7^(f + 0)), 
Csie) = (7r(r)7J(r + 0)), 



(23) 
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where the angled brackets denote the average over all galaxy pairs 
separated by an angle 0. The superscript, r denotes rotated shear 
components which are equivalent to 71 and 72 in a rotated coor- 
dinate frame, defined by the line joining the centroids of the two 
galaxies in question. In terms of the elements of the shear matrix 
(equation (3)), these rotated shear components are 

71 = 71 cos{2(j)) + 72 sin(2</)) 

7^ = -7isin(20) -|-72cos(20), (24) 

where is the angle between the original and rotated coordinate 

frames. Once again, the parity invariance of weak lensing predicts 
that the cross-correlation function 03(6) should be zero. A non- 
zero C3 (6*) is, therefore, an indication of residual systematic effects 
present in the data. The correlation functions, C'i{9) & C2{9) are 
related to the underlying convergence power spectrum via 

Ci{e) = ^CT [Mm + i-iy+'Mm] , (25) 

where J„ (x) are Bessel functions. 



2.2.3 Shear variance in cells 

Finally, we can define the variance, a^{d) of the shear field, mea- 
sured in circular cells of radius 9. In terms of the convergence 
power spectrum, this variance is 

= if ."cr(H^)". 

In the following analysis, we have chosen to meastire the shear vari- 
ance in square cells of side length, 0. However, the corresponding 
variance in circular apertures can be recovered, to a good approx- 
imation, by scaling our measurements by l/\/7r (e.g. Bacon, Re- 
fregier & EUis 2000). 

Having considered the shear field and its statistical properties 
we now turn to our data set. In Sections 3 & 4 we describe the 
COMBO- 17 dataset and our method used in producing the high 
quality deep R-band images from which we measure the shear field. 
In Sections 5 & 6 we shall use this data to determine the statistical 
properties of the observed shear field. 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 The COMBO-17 survey 

The observations analysed in this paper have been undertaken as 
part of the COMBO-17 survey (Wolf et al. 2001). This survey is be- 
ing carried out with the Wide-Field Imager (WFI) at the MPG/ESO 
2.2m telescope on La Silla, Chile. The survey currently consists 
of five 0.5° X 0.5° fields totaling 1.25 square degrees with obser- 
vations taken in five broad-band filters (UBVRI) and 12 narrow- 
band filters ranging from 420 to 914 nm. The chosen filter set facil- 
itates accurate photometric redshift estimation (iTz « 0.05) reliable 
down to an Ji-band magnitude of 24. During the observing runs, to 
facilitate accurate weak lensing studies, the best seeing conditions 
were reserved for obtaining deep _R-band images of the five fields. 
It is these 7?-band images, along with the photometric redshift ta- 
bles, that we make use of in this analysis. 



3.2 Initial data processing 

The initial reduction of the data proceeded along the lines of that 
described in Gray et al. (2002, hereafter GTM+). The WFI instru- 
ment consists of a 4 X 2 array of 2048 x 4096 pixel CCDs. With a 
pixel scale of 0.238" , the resulting total field of view (FOV) of the 
WFT is 0.56° X 0.55°. The standard COMBO-17 pre-processing 
pipeline (Wolf et al. in prep.) produces mosaics of eight 2K x 4K 
chip images which have been debiased, corrected for non-linearity, 
normalized, flat-fielded and cleaned of cosmic rays. The resulting 
mosaics have been simply skewed as an approximate correction 
for rotational mis-alignments between the chips. We have removed 
this skewing of the mosaiced images in order to restore the origi- 
nal chip configurations on the detector The individual chip images 
were then extracted from the mosaics so that astrometric calibra- 
tions could be applied to each of the chip images individually. 



3.3 The precise astrometric solution 

After extraction from the mosaics, catalogues of objects were cre- 
ated for each chip exposure using SExtractor (Bcrtin & Arnouts 
1996). Using pointing information from the image headers, an ap- 
proximate transformation was calculated to convert the SExtractor 
(a;, y) pixel coordinates of the objects to celestial coordinates {a, S) 
for each exposure. These objects were then matched with objects 
from the digitized SuperCOSMOS Sky Survey (SSS)^ (Hambly et 
al. 2001) to within a tolerance of 5". The resulting coordinate pairs 
were then used to iteratively calculate a linear astrometric solution 
for each chip image with the optical axis being used as the tangent 
point for projection. Typically around 300 objects per chip image 
were used for the final astrometric fits which had rms residuals of 
~ 0.2" (c.f. the pixel scale of 0.238"/pixel). Including higher order 
terms and/or a radial distortion term in the astrometric solution has 
been investigated by GTM+. Although, in some cases, the fit was 
improved when a radial distortion term was included, they found no 
improvement for the instrument as a whole. Details concerning the 
artificial shear introduced by such a radial distortion are discussed 
in GTM+ who deduced a radial distortion of 5r/r ~ 0.025% and 
a resulting instrumental shear pattern with amplitude 7 < 0.0001. 
This level of distortion is clearly negligible in comparison to a typ- 
ical ACDM cosmic shear signal of 7 ~ 0.01 and so we have used 
the simple linear solution to perform astrometry on all the chip im- 
ages. Further tests of the artificial shear introduced by both the WFI 
instrument and the telescope optics are described in the next sec- 
tion. Having produced linear astrometric fits for all the chip images 
in each of the five fields, these fits were then used to register the im- 
ages to the same coordinate system. The images were thus aligned 
and combined, using 3a bad pixel rejection with weighting by ex- 
posure time and scaling by the median pixel value. All bad columns 
and pixels were removed during the co-addition procedure due to 
the large number of images combined. The five resulting 8K x 8K 
images were trimmed to remove the under-sampled edges of the 
fields resulting in five 8192 x 8192 pixel (32.5' x 32.5') images. 



t Full access to the SSS is available via the World Wide Web at URL 
http : / /www-wf au . roe . ac . uk 
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Figure 1. Instrumental shear pattern for the 8 chips of the WFI instrument as measured from the observed positions of the same objects on three dithered 
frames. The length of each vector represents the magnitude of the instrumental shear at the indicated position on the chip and the direction indicates the 
orientation of the distortion. For comparison, a 0.1% shear is also shown on the same scale. 



3.4 Artificial shear introduced due to instrument and image 
co-addition 

We have investigated the level of distortions introduced by the tele- 
scope optics and the WFI instrument by comparing the location 
of the same objects on dithered chip exposures. As documented in 
Bacon et al. (2000, hereafter BRE), the respective positions of the 
same object on two dithered frames, / and /' can be expressed as 

x^' -x^ ~^{x^ -x^'), ill) 

where x^ (x^ ) is the position of the object on frame / (/'), x^ 
{x^ ) is the position of the centre of frame / (/') and is the 
distortion matrix, which may be expressed in terms of observables 
as 

^^f^ + ^i V (28) 

Here, ft and 7, are the spurious convergence and shear introduced 
by the geometrical distortions. By measuring the positions of the 
same objects on three dithered frames, the resulting two equations 
of form (27) can be solved for k and 74. We have used this method 
to map the instrumental distortion across the 8 component chips of 
the WFI and the resulting shear pattern is shown in Fig. 1 along 
with an indication of a 0.1% shear signal. We find, in agreement 
with GTM+, that the induced shear due to telescope and instru- 
ment distortions is indeed negligible (7 < 0.2%) over the entire 
field compared to the typical weak lensing signal one would wish 
to measure. 



We have also estimated the tmcertainty introduced in the ellip- 
ticity measurements of objects on the combined images as a result 
of the co-adding procedure. Having registered the original expo- 
sures to the same coordinate frame, we then calculated the disper- 
sion in the positions of many objects from the measured positions 
of those objects on all frames. We found a mean rms of « 0.5 pix- 
els in the positions of objects. Using this mean uncertainty in the 
position of objects, we can simulate the effect of image co-addition 
on a circular source of size ^0.7 arcsec, similar to the inean seeing 
disc size on the individual frames. We simulate the stacking of 78 
frames (c.f. no. of frames for CDFS field - see Table 1) where we 
apply to each image an offset, dx = + 5x where da;*'"'' is 

the true offset position of the object and 5x is a vector whose mag- 
nitude is taken from a Gaussian distribution centred on 0.5 pixels 
and whose orientation is taken from a uniform distribution between 
and 27r. The size and ellipticity of the resulting stacked object 
was then computed using the formalism described in section 5.3 of 
BRE. Fig. 2 shows a histogram of the resulting ellipticity measure- 
ments for 1000 such simulations demonstrating that the shear aris- 
ing from co-addition is again well below the sought-after cosmic 
shear signal. The median ellipticity induced in a circular source as 
measured from these simulations was Cm = 0.0037 and the inter- 
quartile range was 5em = 0.002. These numbers are well within 
our error budget for a precise measurement of cosmic shear. 
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Figure 2. Ellipticity introduced in a circular source due to co-addition of 
individual frames as measured from 1000 simulations of the stacking of 78 
frames. 
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3.5 Point spread function corrections and generation of 
object catalogues 

We have used the imcat software (Kaiser, Squires and Broad- 
hurst, 1995; hereafter KSB) to generate object catalogues for the 
five fields, measure shape parameters for all objects and to correct 
for such effects as the isotropic smearing of the point spread func- 
tion (PSF) of objects by the atmosphere and telescope as well as 
any anisotropic smearing introduced by tracking errors, imperfect 
dither alignments and co-addition of the individual frames. We have 
used the hf indpeaks routine to locate objects, followed by the 
getsky, apphot and getshapes routines to estimate the 
local sky background, measure aperture magnitudes and half-light 
radii, and to calculate shape parameters for all detected objects 
(see e.g. KSB; BRE; GTM-I-, for details of how the KSB procedure 
and the imcat software works). 

This procedure produced an image catalogue for each field, 
containing positions, magnitudes, sizes and shapes for all the ob- 
jects detected. We then applied reasonably conservative cuts on im- 
age size, signal-to-noise (S/N) and object ellipticity (r-g > 1.0, v > 
5, e < 0.5) to remove spurious and/or untrustworthy detections. 
Note that the S/N cut we have used (v > 5) is much less stringent 
than that adopted in BRE who used a S/N cut of f > 15. The mo- 
tivation in BRE for using such a conservative cut was the existence 
in their data of a highly significant anti-correlation between their 
measured mean shear values and the corresponding mean stellar 
ellipticities for a S/N cut of > 5. We have searched for this effect 
in our data on several dififerent smoothing scales and have found 
no evidence of correlations between our shear and stellar ellipticity 
values. For example. Fig. 3 shows the mean shear components in 
cells plotted against the mean uncorrected stellar ellipticity compo- 
nents for a cell size of 8 x 8 arcmin. This is the same cell size as 
that shown in BRE's fig. 7 where a clear anti-correlation is present. 
It is clear from Fig. 3 however, that this effect is absent from our 
data. We have looked for such an effect on scales ranging from 1 to 
30 arcmin and have found none. We have therefore used a z/ > 5 
cut in our analysis. We also note that this is the same S/N cut as 
used in GTM-l- for their supercluster analysis of the A901 field. 

Having applied these cuts to the data, the imcat software 
was then used to correct the measured shapes of the galaxies for 



Figure 3. Mean shear components in cells plotted against the mean stellar 
ellipticity components for a cell size of 8 x 8 arcmin. The upper panel is 
for a S/N threshold of u > 15 and the lower panel is for a threshold of 

u > 5. The anti-correlation of 7^ with e* found in the analysis of BRE is 
not present in the COMBO-17 data for either of the S/N cuts. 

the effects mentioned above. In particular, we correct the galaxy 
shapes for (a) anisotropy and (b) circularization. Details of the cor- 
rection scheme used by imcat have been described elsewhere (e.g. 
KSB; Luppino & Kaiser 1997; Hockstra ct al. 1998; Kaiser et al. 
1999) and will not be repeated here. Note that these corrections 
were done on sections (size 8' x 16') of the fields rather than the 
full 32.5' X 32.5' fields themselves. This was done to ensiu'e that 
an accurate model of the PSF distortions was obtained. Perform- 
ing the corrections on the full fields would possibly lead to residual 
PSF distortions in the dataset due to inaccurate modelling of the 
stellar ellipticities over such a large field. 

A demonstration of the anisotropic PSF correction at work for 
the Sll field is given in Figs. 4 and 5. The smoothly varying stel- 
lar ellipticity pattern apparent in the stars before correction (Fig. 5, 
left-hand side) is successfully removed to produce residual stellar 
ellipticities with essentially random orientations (Fig. 5, right-hand 
side) and mean residuals of I (5ei I ~ 5x 10"'' and |(5e2l ~ 2x10"''. 
Fig. 4 shows the stellar ellipticities after correction randomly dis- 
tributed about ei = 62 = 0. Note however that these plots are for 
illustrative purposes only. Rigorous tests of the systematics intro- 
duced by residual spurious ellipticities are performed later in Sec- 
tion 5, where we measure the star-galaxy cross-correlation fimc- 
tions for the dataset. 



4 THE COMBO-17 FIELDS 
4.1 Content of the fields 

The five fields observed are quite different from one another in 
terms of their content. Three are blank fields; one (CDFS) was cho- 
sen so it would overlap with the Chandra Deep Field South; another 
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Figure 5. Stellar ellipticity pattern across the Sll field before {left) and after {right) correction for PSF anisotropy. A 10% shear signal is also shown for 
comparison. The horizontal and vertical lines show the sections used for the stellar ellipticity polynomial model fitting. 



Table 1. Properties of the five COMBO-17 shear catalogues. In the sixth column we quote the median R-band magnitude of each catalogue. In the case of the 
SGP and FDF fields, the median mpj^ listed is an approximate value as magnitudes for these fields have been roughly calibrated using the APM galaxy survey. 
The mean galaxy density values (rigai) are calculated as rigai = Nf/Af where Nf and Af are the total number of galaxies in the final shear catalogues and 
the total useable area of the fields respectively, after excluding regions contaminated by bright stars, diffraction spikes and ghosting (see Fig. 6). The FWHM 
values Usted are the average FWHM of stars on each field as measured from the final co-added images. The final column indicates whether photometric redshift 
information is currently available for the field. 



Field RA(2000) Dec.{2000) R-band Exposures Total Exp. Time Median mR ngai (arcmin^^) FWHM Redshift info.? 



CDFS 


03"^ 32™ 25= 


-27° 


48' 


50" 


42 X 500s -f 36 X 420s 


36120s 


24.4 


37.5 


0.81" 


Yes 


SGP 


OQh 45™ 56= 


-29° 


35' 


15" 


42 X 500s 


21000s 


23.6 


36.2 


0.81" 


Pendin 


FDF 


Oih 05m 49s 


-25° 


51' 


42" 


20 X 500s -1- 5 X 400s 


12000s 


23.7 


28.4 


0.82" 


No 


Sll 


111" 42™ 58= 


-01° 


42' 


50" 


39 X 400s -1- 5 X 500s 


18100s 


23.7 


27.7 


0.78" 


Yes 


A90I 


09^ 56™ 17= 


-10° 


01' 


25" 


36 X 500s -1- 8 X 600s 


22800s 


24.0 


36.1 


0.76" 


Yes 



(SGP) is centred on the South Galactic Pole, and the third blank 
field is the shallower FORS deep field (FDF). A fourth, randomly 
selected field (Sll) contains a fairly large cluster (Abell 1364) at 
a redshift of z = 0.11. The fifth field (A901) was chosen so as to 
include a supercluster system (Abell 901/902) at redshift z = 0.16 
and the COMBO-17 observations of this field have already been 
the subject of an extensive weak lensing analysis, details of which 
are given in GTM-i-. Details of the R-haad observations for the five 
fields are given in Table 1. Note that we do not currently have the 
full 17-band photometric information for the SGP or FDF fields. 

The significant mass concentrations present in the two fields 
containing clusters (Sll & A901) raise the question as to whether 
these fields should be included in a cosmic shear analysis at all. The 
significant difference between these fields is the pre-selection of 
A901 to include a supercluster system. Strictly speaking, therefore, 
the A901 field is not a randomly selected piece of sky and is un- 
likely to be a fair representation of the Universe. We have therefore 
chosen not to include this field in our final parameter estimation 
and shear power spectrum calculations. We have also experimented 
with including/excluding both the Sll and A901 fields from the 



analysis. For both the shear power spectrum reconstruction and the 
parameter estimation, we find that the A901 field biases our results 
significantly, whereas the Sll field is consistent with the dataset 
as a whole (see Section 6.3). The cell-averaged shear distributions 
for the five fields are shown in Fig. 6, showing regions in the fields 
which have been masked out due to contamination by bright stars, 
diffractions spikes and ghosting. 



4.2 Including redshift information 

As mentioned earlier, the COMBO-17 survey has the unique advan- 
tage, in terms of weak lensing surveys, of having accurate photo- 
metric redshift measurements for ~ 90% of objects detected down 
to a magnitude of uir < 24 (for a detailed description of the pho- 
tometric methods used to assign redshifts to the galaxies, see Wolf 
et al. 2003). For these galaxies, the estimated typical uncertainty 
in redshift is az 0.05. At the time of writing, photometric data 
is available for three of our five fields (CDFS, A901 & Sll). We 
have used this information to estimate, directly from the data, the 
redshift distribution and in particular, the median redshift, Zm of 
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(arcmin) 
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FDF 

X - 7 i \ 



10 20 

(orcmin) 



A901 



10 30 

(arcmin) 




-'■-'-JX' - \ -7/- 
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10 30 

9 (arcmin) 



10% shear 



Figure 6. The shear distribution on the five COMBO fields. The fields are (clockwise from top left) CDFS, SGP. FDF, SI 1 and A901. The shear measurements 
for each field have been binned into 30 X 30 pixels, giving a pixel scale of ^ 1 arcmin. For each pixel, the length of the vector drawn is proportional to 
the magnitude of the mean shear in that pixel and its direction indicates the orientation of the shear. A 10% shear signal is also shown for comparison. The 
apparent holes in the distributions are regions which have been masked out due to bright stars, diffraction spikes etc. 



0.10 



0.05 - 
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e, 

Figure 4. Ellipticity distribution of stars in the Sll field before (grey 
crosses) and after (black crosses) correction for PSF anisotropy. Mean 
stellar ellipticity components of eT = —0.025 ± 0.013 and 'e-i = 
—0.007 ± 0.011 are removed to leave an essentially randomly orientated 
stellar ellipticity distribution with mean residual components of &e\ = 
-0.0005 ± 0.010 and 5el = -0.0002 ± 0.009. 



the lensed source galaxies in the survey. For the Zm calculation, we 
have only included those galaxies used in the final weak lensing 
analysis. This subset of galaxies does, in fact, include galaxies with 
magnitudes mn > 24 for which the redshift estimations are unreli- 
able. To account for this, we must extrapolate the redshift distribu- 



Table 2. Measured median magnitudes and redshifts (Zm) for the three 
fields, CDFS, A901 & Sll for different limiting magnitude cut-offs. Also 
listed are the number of galaxies (A^g^i) used in the median calculations 
and the completeness of each magnitude-limited sample (i.e. the fraction of 
galaxies in the sample for which accurate redshifts have been obtained). Ex- 
trapolating the median rrtjs - Zm relation (see Fig. 7) to the median magni- 
tude of the combined shear catalogues for all fields (= 24.0), we infer a me- 
dian redshift for the COMBO-17 weak lensing survey of = 0.85it0.05. 



Limiting rriR 




Median 




Completeness 


< 20.0 


1485 


19.15 


0.184 ±0.004 


87% 


< 21.0 


3267 


20.13 


0.285 ± 0.030 


87% 


< 22.0 


7316 


21.14 


0.404 ± 0.026 


88% 


< 23.0 


16387 


22.14 


0.520 ±0.023 


87% 


< 24.0 


31649 


22.96 


0.673 ± 0.092 


76% 



tion beyond m_R = 24. For the three fields for which photometric 
data is available, accurate redshifts have been found for 31649 out 
of the 83514 galaxies used in the final shear analysis (~ 38%). 

For a number of different limiting magnitudes, we have mea- 
sured the median magnitude of the galaxies and the correspond- 
ing median redshift {zm) of those same galaxies. Details of these 
measurements are given in Table 2. We take our median redshift 
measurements as lower limits for the true median redshift as the 
magnitude-limited samples are not complete in redshift (see Table 
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2) and those galaxies without an assigned redshift are most likely to 
be at a higher redshift than our measured median value. To calcu- 
late an upper limit for our Zm estimates, for each magnitude limited 
sample, we have placed those galaxies without redshift measure- 
ments at z — oo and have re-calculatcd z,-,-, . We take our final me- 
dian redshift estimate to be simply the midpoint between our upper 
and lower limits. 

We have then extrapolated the median relation to 

find an estimate of the median redshift of our galaxy sample as 
a whole, which has a measured median magnitude of 24.0. The 
median niR — z,n relation is shown in Fig. 7. Here we plot both the 
COMBO- 17 data points and data from the HDF redshift survey of 
Cohen et al. 2000, with which we are fully compatible. The curve 
plotted is the best-fit quadratic model to the COMBO- 17 data and 
is given by 



Zm = 2.53 — 0.33 mR,m + 0.01 ttir , 



(29) 



where mn,Yn is the median R band magnitude of the galaxy sample. 
From this, we estimate a median redshift of Zm = 0.85 ± 0.05 for 
the COMBO- 17 weak lensing survey. 

There is, of course, more information in the measured redshift 
distribution than just the median value. The advantage of including 
photometric information in a weak lensing analysis lies in reduc- 
ing or eliminating uncertainties in the redshift distribution of the 
source galaxies when it comes to comparing weak lensing mea- 
surements with those predicted from theory (see Sections 6 & 7). 
To this end, we have included the measured redshift distribution for 
our comparison with theoretical models. There are, of course, many 
more galaxies for which we have no redshift estimate than there are 
galaxies with reliable measurements. This situation will improve as 
the COMBO- 17 survey nears completion, but for the purposes of 
this analysis and predicting the expected weak lensing signal for 
different cosmological scenarios, we need to assign redshifts to the 
galaxies with unknown redshifts. We do this by distributing these 
galaxies in z according to (Baugh & Efstathiou 1991) 



dN 
dz 



pz^ 



exp 



(i)' 



(30) 



with = 1.5 and Zi, = 2;m/l-412. A distribution of the form, 
equation (30) has a median redshift of Zm- We have tuned the value 
of Zm in equation (30) to ensure that the final total redshift distri- 
bution has a measured median value of Zm = 0.85. 

After assigning values to the galaxies with unknown redshifts 
according to equation (30), we now have the final n{z) distribution 
which we can use for making predictions for the weak lensing sig- 
nal expected in the COMBO-17 survey. This final n{z) distribution 
is shown in Fig. 8. 



X Combo- 17 

A HDF (Cohen et al. 2000) 




20 



22 



24 



median m,. 



Figure 7. The dependence of the median redshift (zm) of the COMBO-17 
galaxies on the median magnitude of the galaxies as measured from the 3 
fields for which we have photometric data. A simple quadratic model (solid 
curve) for the — Zm relation is also plotted for comparison. For our 
final weak lensing dataset, which has a median R-band magnitude of = 
24.0, we infer from this trend a median redshift of Zm = 0.85 ± 0.05. 




5 CORRELATION FUNCTIONS AND SHEAR 
VARIANCE ESTIMATORS 

It is useful to estimate the correlation functions of the shear field 
both as a first approach at measuring the strength of the cosmic 
shear signal, but also as a test for unwanted systematics in the data. 
These tests have now become a standard part of shear analysis (e.g. 
Van Waerbeke et al. 2001; Pen, Van Waerbeke & Mellier 2002; 
Bacon et al. 2002). Measuring correlation functions also allows a 
direct comparison of our measurements with previous cosmic shear 
studies. Here we present both an unweighted correlation function 
analysis and a minimum variance weighted shear variance analy- 
sis. Note that, for the remainder of the paper, we have excluded 



Figure 8. The final redshift distribution (solid line) used for making model 
predictions to compare with the observed weak lensing signal from the 
COMBO-17 dataset. This composite distribution has a median redshift of 
Zm = 0.85 and is composed of a measured n{z) distribution (from the 
three fields, CDFS, Sll and A901; shown here as the dotted line) and a 
parameterised model (equation (30); dashed line) for the galaxies with un- 
known redshifts. 



the A901/2 supercluster field from our measurements except where 
specifically indicated. 
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5.1 Correlation analysis 



To estimate the unweighted correlation functions, we follow Ba- 
con et al. (2002). For the purposes of the correlation analysis, we 
have divided each of the COMBO- 17 fields into eight, chip-sized 
sections. This allows for a better estimate of the field-to-field co- 
variance between our measurements than would be possible using 
the original 0.5° x 0.5° fields. The correlation functions for the in- 
dividual sections can be measured by averaging over galaxy pairs: 

Us ^-^ 

i=l 

C2AO) = - V75(fO75(n + 0), 

Us ^-^ 

i=l 

CzAQ) = —yYl{r^hKr^ + 0), (31) 

i=l 

where is the number of pairs of galaxies in each section and 
where we average over all galaxy pairs separated by an angle, Q 
(c.f. equation (23)). The signal averaged over all sections is simply 



1 



(32) 



where Ns is the total number of sections. The covariance of the 
correlation function measurements is given by 



cov[C.(6»)Cj(6»')] 



(33) 



Fig. 9 shows the correlation functions, Ci{9) and C2{9) af- 
ter averaging over all sections of the four fields, CDFS, Sll, SOP 
and FDF. As can be seen, our data agrees well with the expected 
correlations, for both Ci and C2, for a — 0.3 ACDM model 
with a power spectrum normalisation of as — 0.8 (see Section 7 
for our measurements of fi,„ & as from COMBO- 17). These the- 
oretical curves have been calculated using equation (25) where we 
have used the halofit model of Smith et al. (2002)^ to calculate the 
non-linear matter power spectrum. We have also input our com- 
bined redshift distribution (Fig. 8) into the calculation. We should 
point out that on smaller scales than those shown on Fig. 9, i.e. 
6 ^ 0.3 arcmin, we have measured an inconsistently large sig- 
nal (^ 5 X 10"*) for the Ci correlation function. We suspect that 
this greatly enhanced signal is probably due to either systematic 
effects on very small scales due to imperfections in the correction 
for PSF anisotropy, or alternatively, the increased signal may be 
caused by intrinsic alignments dominating at such small angular 
separations. Either way, we exclude this single data point for our 
parameter estimations in Section 7. We also note that the errorbars 
on Ci{6) and 02(6) may be slightly underestimated on the very 
largest scales due to correlations between neighbouring sections of 
the same field. Fig. 10 shows the cross-correlation between shear 
components, C-i{9) which is consistent with zero, as expected. 

In Fig. 11, we show the star-galaxy cross-correlations for 
our data. The systematically-induced cross-correlation between the 



■t- the halofit code is publicly available at URL 

http : / /www . asl . chem . nott ingham .ac.uk/~res/ software 



X 
CO 



X 



n =0.3, iT„=0.8, z =0.85 

m o m 

« Cj(6), prediction: 

» Cg(e), prediction: 




o 



G (arcmin) 

Figure 9. The unweighted sheai' congelation functions Ci{6) and C2{9) 
from the COMBO- 17 data. The solid dark line is the expected Ci corre- 
lation for a Qrn. = 0.3 ACDM model, normalised to erg = 0.8, while 
the lighter dashed line is the expected C2 signal. The C2 points have been 
slightly displaced horizontally for clarity. Note that we measure an incon- 
sistently large Ci signal (> 5 X 10~*) on the very smallest scales (9 < 0.3 
arcmin. not shown) which, we suspect is due to residual systematic effects 
and/or intrinsic alignments. 



galaxies (corrected for PSF anisotropy) and the uncorrected stars is 
defined by (Bacon et al. 2002): 



(34) 



where i = 1, 2. A correlation here would indicate a problem in 
the correction of systematic image distortions. Again we find that 
no correlations are detected at scales 9 ^ 0.3 arcmin, indicating 
that our PSF correction procedure has worked well. Note, however 
that the star-galaxy cross-correlations are significant at scales 9 < 
0.3 arcmin, indicating a possible systematic problem at these very 
small scales - this motivates our exclusion of the smallest-scale 
Ci {9) data point for our parameter estimations in Section 7. 

We can use our correlation function measurements to com- 
pare our results with those of previous cosmic shear studies. We 
perform such a comparison by looking at the total correlation func- 
tion, C{9) = Ci{9) + C2{9) for each data set. In Fig. 12, we 
plot C{9) as measured from COMBO- 17, along with the most re- 
cent results from four other groups (Van Waerbeke et al. 2001; Ba- 
con et al. 2002; Refregier et al. 2002; Hoekstra et al. 2002). Note 
that we have scaled each groups' results to a median redshift of 
Zs = 0.85 using the scaling suggested by the numerical simula- 
tions of Barber (2002), i.e. C{9) oc z^^. We have also, in the case 
of Refregier et al., and Hoekstra et al., assumed that a?^{9) ~ C{9), 
where a?^{9) is the shear variance statistic measured by these two 
groups. That this is a reasonable assumption can be seen by com- 
paring the as — 0.9 models in Figs. 12 and 13. Fig. 12 shows that, 
beyond a scale of 6* ~ 1 arcmin, all groups are in broad agreement 
with each other, although the COMBO- 17 signal does seem to be 
slightly lower in amplitude. At scales smaller than 1 arcmin, com- 
htmJparisons are more difficult to make because of the larger error bars 
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Figure 10. The cross-coiTelation of the shear components 03(6), as mea- 
sured from the COMBO-17 dataset. The signal is consistent with zero on 
all scales in agreement with theoretical predictions. 



9 (arcmin) 

Figure 12. The total correlation function, C{8), as measured from 
COMBO-17, along with the most recent cosmic shear measurements from 
the four other groups indicated. Beyond a scale of 1 arcmin, the measure- 
ments are in broad agreement. The correlation function predictions for a 
flat ACDM cosmology, for three values of the power spectrum normalisa- 
tion (from top to bottom; erg = 1.0, 0.9, 0.8) are also plotted. 
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Figure 11. The star-galaxy cross-correlation functions (equation (34)) for 
the COMBO-17 data. The Ci and C'2 points have been slightly offset hor- 
izontally for clarity. These measurements are a strong test of how success- 
fully the galaxies have been corrected for distortions introduced by PSF 
anisotropy. The residual coiTelations between galaxies and stars are con- 
sistent with zero on all scales. In the worst-case scenario, represented by 
the upper end of the error bars shown, our correlation function measure- 
ments would include < 10% contamination from residual distortions left 
over from the PSF corrections applied. 
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Figure 13. The minimum variance shear variance as measured from the 
COMBO-17 dataset. Also plotted are predictions for the shear variance for 
a Qm = 0.3 ACDM cosmology for three different normalisations of the 
matter power spectrum. 



involved. Indeed, it is the measurements at the larger scales (0^1 
arcmin) that provide our best constraints on cosmological models. 
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5.2 Shear variance 

We have also applied the minimum variance estimator of Brown et 
al. (2002) to measure the shear variance in square cells in excess of 
the noise: 



.2 ,n\ _ Zlcell Wcell(e^ell - ^cell) 



Eoell "'crfl 



(35) 



where Cceii is the cell averaged elUpticity and Wceii is a minimum 
variance weight; 

1 



Wcell ■ 



(36) 



Here TVceii is the noise in the cell measured from the data due to the 
intrinsic dispersion in galaxy ellipticities, 



n 

JVceiiW = ^5^(e?(fi)+e2(ri)), 



(37) 



where n is the number of galaxies in a cell and is the pre- 
dicted shear variance for square cells which we calculate for a 
ACDM model using x equation (26). Varying the cell side 
length yields the cell variance as a function of scale, aien3(^')- The 
error in equation (35) is given by 



Var[a-iens] = 



(38) 



In Fig. 13, we plot our shear variance measurements along 
with the predicted shear variance signal for a ACDM model, 
and for three values of the power spectrum normalisation, as = 
0.8, 0.9 & 1.0. Immediately, one sees that the shear variance statis- 
tic is measuring a somewhat higher signal than the correlation func- 
tion on scales < 2 arcmin. We suspect that this is the same effect as 
that seen in the correlation functions at very small scales, i.e. that at 
these scales residual systematics and/or intrinsic galaxy alignments 
become important. Very recently, methods have been developed 
for separating the intrinsic and lensing signals (Heymans & Heav- 
ens 2003; King & Schneider 2002) by effectively down- weighting 
physically close pairs of galaxies when calculating the correlation 
functions. Such methods are dependent on accurate photometric 
information being available for the individual galaxies. In this re- 
spect, the COMBO-17 survey is ideal for the separation of the in- 
trinsic and lensing signals, and we are currently applying such an 
analysis to the survey. We note, for the purposes of the current lens- 
ing £inalysis, however, that intrinsic alignments are unlikely to con- 
tribute significantly to our signal beyond scales of a few arcmin in 
the shear variance or beyond ~ 1 arcmin in the correlation func- 
tions. Indeed, the few observational studies of intrinsic alignments 
carried out to date (Pen, Lee & Seljak 2000; Brown et al. 2002; 
Lee & Pen 2002) support arguments from both theoretical consid- 
erations (Crittenden et al. 2001; Catelan et al. 2001; Mackey et al. 
2002) and numerical simulations (Heavens, Refregier & Heymans 
2000; Croft & Metzler 2001) that predict the intrinsic contribution 
to be no more than ~ 10% for deep surveys such as COMBO-17. 
We do, however, urge caution when interpreting these statistics on 
scales ^ 1 arcmin. 



6 COSMIC SHEAR LIKELIHOOD ANALYSIS 

Having tested our data for sources of systematic errors and mea- 
suring the cosmic shear signal by standard methods, we now apply 
a new maximum likelihood analysis to estimate the cosmic shear 



power spectrum from the COMBO-17 dataset. This section con- 
tains the key results of our paper 

6.1 Likelihood procedure 

Our approach is based on the prescription of Hu & White (2001, 
hereafter HW), who use a maximum hkelihood method to recon- 
struct the three power spectra, C^", C'l'^ and ^, directly from 
pixelised shear data. More precisely HW proposed reconstructing 
Cg'^ as a series of step-wise "band powers", extracted from the 
data via an iterated quadratic estimator of the maximum likelihood 
solution. This approach, which is similar to methods used in anal- 
yses of CMB polarisation fields (e.g. Tegmark & de-Oliviera Costa 
2002), has the advantage that it automatically accounts for irregu- 
lar survey geometries and produces error estimates which include 
sampling variance and shot noise. In addition this approach can ac- 
count for the effects of pixelisation. 

HW have tested the maximum likelihood estimator on both 
Gaussian realisations of a ACDM power spectrum and on N-body 
simulations. In both cases, the estimator performs well, recovering 
the input power accurately with error estimates from the Fisher ma- 
trix showing excellent agreement with run-to-run errors. We also 
investigate the method with our own simulations in the next sec- 
tion. 

Another advantage of this approach for reconstructing the 
shear power spectrum is that it provides a simple method for per- 
forming a decomposition of the signal into curl ((3) and curl-free 
(k) modes. The weak lensing shear power spectrum is predicted to 
be completely curl-free in the absence of significant lensing from 
gravitational waves (Stebbins 1997; Kamionkowski et al. 1998) al- 
though the predictions for intrinsic galaxy aUgnments are less cer- 
tain (e.g. Crittenden et al. 2001, Mackey et al. 2002). The Klf5 de- 
composition, therefore, represents a useful method for the detec- 
tion of non-lensing artefacts (e.g. intrinsic alignments, systematic 
effects) in the data. 

If we write our data as a vector, 

d = (7i(T'i),72(ri),--- ,7l(^'n),72(rn)), 
then the likelihood function is 

-21nL(C|d) = d*C"^d-|-Tr[lnC], (39) 
where 



{dd') 



(40) 



is the data covariance matrix, and we assume uniform priors. We 
can interpret the data covariance matrix as the sum of the shear 
covariance matrix (equation (21)) and a noise term 



TV — J;^T 



(41) 



which we measure directly from the data. Here, / is the identity 
matrix, 7ri„s is the intrinsic dispersion of galaxy ellipticities within 
a pixel and iVpix is the pixel occupation number 

Following HW, we maximise the UkeUhood as a function of 
the model parameters. Here, our model parameters are just the band 
powers of the three power spectra, C/'"' , Cf and C^'' and we per- 
form the maximization iteratively with a Newton-Raphson scheme. 
That is, from an initial guess of the band powers, 6i, a new estimate. 



'/ = ei + AOi 



(42) 



is made for the band powers where we have adjusted our previous 
estimate by 
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AOj = 2 



{ de, ) 



dlnL\ I \nL 



(43) 



Here, we can replace the second derivative of tiie likeliiiood by its 
expectation value, 



' InL 



(44) 



and at the same time, in the limit where the likelihood is sufficiently 
Gaussian in the parameters, we can use the Fisher Information ma- 
trix (e.g. Tegmark, Taylor & Heavens 1997), 



(45) 



to estimate the uncertainties on the band power measurements. That 
is, we approximate the covariance matrix of our band power esti- 
mates as 



(46) 




1000 



6.2 Testing the likelihood on simulations 

Although HW test the likelihood reconstruction of the various 
power spectra on both Gaussian and N-body simulations, we have 
also conducted our own simulations of the likelihood reconstruc- 
tion. We do this because of the small size of our data fields (30 x 30 
arcmin) relative to HW's simulated fields, which are about 50 times 
larger in area. 

We have applied the maximum likelihood estimator to one 
hundred 30 x 30 arcmin fields where the pixelised shear distri- 
bution in each field is a Gaussian realisation of a shear power spec- 
trum calculated for a ACDM cosmology and for all source galaxies 
lying atz — 1. Following HW's approach, we have added Gaussian 
distributed noise to each of the 20 x 20 pixels (pixel size, 1.5 x 1.5 
arcmin) of our shear distribution according to equation (41) where 
we have taken 7rms = 0.4 for the intrinsic dispersion in galaxy 
ellipticity components. For the number of galaxies in each cell, we 
employ a mean galaxy density of rig = 32.0 arcmin"^, similar to 
the measured galaxy density of our dataset (see Table 1), and we 
calculate the pixel occupation number with A^'pix = rigB^^^ where 
Spix is the pixel side length in arcmin. The results from all the sim- 
ulations are summarised in Fig. 14 where we plot the recovered 
band powers of the three power spectra, Cf*, C^^ and C^^, aver- 
aged over all simulation runs. Here, we plot both the run-to-run er- 
rors and the errors as estimated from the Fisher matrix which agree 
well with one another, supporting the use of the Fisher matrix to es- 
timate errors on the band power measurements. We see from these 
simulation results that the maximum likelihood reconstruction re- 
covers the input power very accurately for Fourier modes in the 
range 1000 ^ £ Si 5000. There is some indication from the sim- 
ulations, however, that the method may slightly over-estimate the 
power for i?— modes outside this range. However, this discrepancy 
is much smaller that the precision to which we measure the signal. 
Finally, we note that the likelihood reconstruction estimates of the 
/3-mode power spectrum, and the cross-correlation, C^^ are 
both consistent with zero on all scales, as expected for the input 
Gaussian shear model used to create the shear distributions. 



Figure 14. The recovered band powers of the three spectra, C^", C^^ and 
C'^^ as averaged over 100 Gaussian realisations. Error bars are shown as 
calculated from the run-to-run variations (smaller terminal ends) and also 
as estimated from the Fisher matrix (larger terminal ends). The solid curve 
is the input ACDM power spectrum used to create the simulated shear dis- 



tributions. The C^^ and C^*^ points have been slightly offset horizontally 
for clarity. 



6.3 Maximum likelihood results 

6. 3. 1 Shear power spectra 

Having tested the method on simulations, we now apply the maxi- 
mum likelihood reconstruction to our five data fields, CDFS, SGP, 
FDF, Sll and A901/2. As these fields are widely separated on the 
sky we can treat them as independent and optimally combine them 
afterwards. 

Fig. 15 shows the results of estimating Q'"', (7"'' and Cf 
for the five fields, in 5 band powers. Also plotted is a model curve 
for the ACDM model with f^A = 0.7, fi™ = 0.3, h = 0.68, 
normalised to erg = 0.8 (see Section 7). The linear power spectrum 
has been transformed to a nonlinear one using the halofit formulae 
of Smith et al. (2002) and, again, we have included the composite 
redshift distribution shown in Fig. 8 in this calculation. 

In general we find that the fields containing clusters, Sll and 
A901/2, both yield higher results than the other fields. However, 
as mentioned already, the results for the Sll field seem consistent 
with the dataset as a whole while those for the A901/2 field clearly 
are not. That the shape is broadly the same for the A901 field is 
to be expected since this part of the nonlinear power spectrum is 
dominated by massive clusters (Cooray & Hu 2001). At the median 
redshift of the sample, the peak of the nonlinear convergence power 
spectrum corresponds to around 3/i^^Mpc, the scale on which we 
expect collapsed clusters to have virialised. 

The other fields, CDFS, the SGP and FDF, are fields without 
large structures and have correspondingly lower amplitude spectra. 
Interestingly on the smallest scales, I ~ 4000, the power is roughly 
the same in all five fields. 

Fig. 16 shows the power from an optimal inverse weighting 
combination of the fields excluding A901/2. Details of the opti- 
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Figure 15. The cosmic shear power spectra estimated from each of the five individual fields in the COMBO-17 dataset. The spectra are for (clockwise, from 
top left) CDFS, SGP, FDF, Sll and A901. Note the much higher power recovered from the A901 supercluster field. The C^^ and C'^^ points have been 
slightly offset horizontally for clarity. 



mally combined power spectra measurements are also given in Ta- 
ble 3. Three of our five band power measurements agree within 
their error bars with the = 0.8 normalised ACDM model plot- 
ted. The high level of power measured at £ ~ 4000 seems to be 
present in each of the fields we have analysed. The largest discrep- 
ancy between the model and the measurements occurs at £ = 2000 
where less than half the power is found. 

In order to interpret the scales that we are probing with our 
power spectrum measurements, we note that, in a flat cosmology, 
the angular size will scale roughly as 

2'kDa 



R 



(47) 



where, for a flat ACDM Universe, the angular diameter distance 
(Da) is well approximated by (Broadhurst, Taylor & Peacock 
1995) 



Da{z) 



(48) 



Ho (l + z)(l + 3/4!^„z)' 

An estimate for the redshift of typical deflectors in the survey is 
Zd ~ Q-^Zm and we have estimated Zm ~ 0.85 (see Section 4.1) 
for the median redshift of our survey. Combining these relations 
(with Q.m = 0.3) gives the following approximate conversion be- 
tween the physical scale being probed and the Fourier variable on 
the sky: 

5130 



R 



(49) 



This relation gives the approximate physical scale indicated at the 
top of Fig. 16. 

We can check our results for systematic effects by estimating 
the power in the /?-/? modes, as well as the cross correlation. 



Table 3. Details of the maximum likelihood reconstructed band powers as 
obtained from the optimal combination of the four fields, CDFS, Sll, SGP 
and FDF. The first two columns list the i'-range of the five band powers. Also 
listed are the detections and uncertainties in each band where we quote the 
band powers. Pa = £{£ + l)C^'/2''r for i = k, (3. The first band only has 
Pkk measurements as the k and /3 powers are essentially indistinguishable 
at these large scales. 



f ■ 


■^max 


Pkk(x105) 


P^^(xl05) 


Pk/3(x105) 


247 


594 


2.92 ± 2.05 






595 


1790 


4.95 ± 2.43 


-1.45 ± 1.26 


2.21 ± 1.45 


1791 


2986 


2.20 ±4.75 


-4.38 ± 3.76 


-1.99 ± 3.27 


2987 


6324 


26.62 ±8.12 


13.85 ±7.34 


0.49 ±5.54 


6325 


13906 


11.57 ±22.39 


7.87 ±22.03 


14.83 ± 14.64 



In all of our band powers the /3-/3 correlation is below the detected 
signal in k and is consistent with zero in all but one band power 
at £ ~ 4000. Similarly the cross correlation is well below our 
measurement of shear power, and is consistent with zero except at 
£ ~ 1000 where a significant detection appears. We conclude from 
the minimal power found in these spectra that our results are not 
strongly contaminated by systematic effects. 
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Figure 16. The cosmic sheai' power spectra from COMBO-17, after ex- 
cluding the A901 supercluster field. Plotted on a linear-log scale are C'^'^ 
(circles), C^^ (crosses) and Cg^ (triangles) in 5 band-averaged band pow- 
ers, as a function of multipole, estimated from the optimal combination 
of a maximum likelihood analysis of the four COMBO-17 fields, CDFS, 
SOP, FDF, and Sll. The error bars are estimated from the Fisher matrix. 
The solid curve is the shear power spectrum expected for a erg = 0.8 nor- 
malised ACDM model. Once again, the C^^ and C^^ points have been 
slightly horizontally displaced for clarity. 
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Figure 17. The correlation matrix (see text for details) of the optimally com- 
bined band power measurements recovered from the maximum likelihood 
reconstruction of C^'^ . The area of each circle is proportional to the degree 
of correlation between bands i and j . Filled circles denote that the bands 
ai'e correlated whereas unfilled circles denote an anti-coirelation between 
the bands. The bands are the same as those plotted in Fig. 16, numbered 1 
to 5, in order from left to right. Note the small values of the off-diagonal 
elements, indicating small (anti-)correlations between different band power 
measurements. 



6.3.2 Covariance matrix of band powers 

In addition to measuring the amplitude of the shear power spec- 
trum, it is also important to consider the correlations between band 
powers. We can quantify how much the bands are correlated with 
one another with the correlation matrix, defined by 



Cori 



CoVi 



(50) 



where Covij is the covariance matrix of the band powers, which 
we measure directly from the data. In Fig. 17, we plot the cor- 
relation matrix of the optimally combined band power measure- 
ments shown in Fig. 16. It is clear from this figure that our band 
power measurements show very little correlation with one another, 
the biggest effect being the slight anti-correlation of neighbouring 
bands. Thus, our maximum likelihood band powers are almost in- 
dependent measures of power. 



6.4 Integral power spectra approximations 

As well as a maximum likelihood approach, more direct methods 
have been proposed for recovering the shear power spectra. We 
have applied one of these direct methods, proposed by Schneider 
et al. (2002, hereafter SvWKM), to our data. SvWKM proposed 
reconstructing C^" directly from a correlation function analysis of 
the data, via an inversion of equation (25). Pen et al. (2002) have 
also suggested a similar approach, but with a somewhat different 
implementation, which they apply to the DESCART dataset. 

The correlation function estimator, which can also be formu- 
lated in terms of band powers, is not dependent on the spatial 



distribution of the shear and, as pointed out by SvWKM should 
not suffer from bleeding of power between bands due to pixelisa- 
tion of the data. SvWKM test their estimator by reconstructing the 
power spectrum from a fiducial weak lensing survey of area, A=25 
deg^ from which they assume they have measured the correlation 
functions for Q" < 9 < 2° . For such a survey, the estimator re- 
covers the input power very accurately for ^-modes in the range 
200 <, 6 <, 2 X 10^ and SvWKM also measure minimal covari- 
ance between their band powers over this range. Over the range, 
2n/8ma,x ^ I ^ 2-K/Omin, the correlation function estimator 
is practically unbiased where 6'min(^max) is the smallest (largest) 
scale at which the correlation functions have been measured. Here 
we apply a modified version of the SvWKM statistic to our dataset. 
If we define two new correlation functions, ^+(^) = Ci{O)+C2{0) 
and ^_ (6) = Ci (0) — C'2 (0), the correlation function estimator of 
SvWKM can be written as 



Ce = 2n 



dee[Kii+{e)jQ{ie) + K2i-{e)Ji{ie)] , (si) 



where Jo,4,{l9) are the usual Bessel functions and Ki & K2 de- 
scribe the relative contribution from the correlators to the 
integral. In order to achieve an optimal combination of 5+ and 
^_ contributions to the integral, SvWKM have constructed a func- 
tion, K\{t) and imposed the constraint, K2 = 1 — K\{tj. How- 
ever, in order to decompose our signal into k and /3 modes, we 
have chosen K\ and K2 to be constants, independent of I. Setting 
K\ = K2 = 1/2 yields a purely K-mode estimator, whereas if 
we set Ki = 1/2 and K2 = —1/2, the estimator should recover 
only the /3-mode power. Following SvWKM, we have formulated 
the correlation function estimator in terms of band powers. 
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Figure 18. The shear power spectrum, C^" and the /3-mode spectrum, Cjf 
as estimated via a modified version (see text for details) of SvWKM's cor- 
relation function estimator The C^^ points have been laterally displaced 
for clarity. Again, we have excluded the A901 supercluster field for these 
measurements. Also shown is the shear power spectrum expected for a 
as = 0.8 normaUsed ACDM model. 



Fig. 18 shows the result of applying the correlation function 
estimator to COMBO-17. Again, we have excluded the A901 su- 
percluster field from this calculation. We note that only the third 
band power on this plot is within the stated range of validity of the 
SvWKM estimator. Indeed, this data point is in good agreement 
with the Q.m = 0.3 ACDM model plotted and the power in /3 for 
this band power is consistent with zero. On both larger and smaller 
scales however, as predicted by SvWKM, the estimates becomes 
unreliable, measuring equal n and (3 modes on large scales. Ide- 
ally, of course, one would measure the correlation functions out to 
as large a scale as possible to extend the range of applicability of 
the estimator. However, for our dataset, we only have correlation 
function measurements in the range 20" ^ 6 12' giving us a 
valid £-range of 1800 ^ t ^ 5000. We conclude, therefore, that 
although reconstructing the shear power spectrum from a correla- 
tion function analysis can be useful for datasets with larger fields, 
for our purposes, the full maximum likelihood analysis is clearly 
the most suitable approach. 



7 COSMOLOGICAL PARAMETER ESTIMATION 
7.1 ^Ira and from COMBO-17 

Having measured various statistics from the data, we are now in a 
position to use both our measured shear correlation functions and 
our reconstructed shear power spectrum estimates to obtain a joint 
measurement of the normalisation of the mass power spectrum (Tg, 
and the matter density Q,m ■ We can achieve this by fitting theoret- 
ical shear correlation functions and power spectra, calculated for 
particular values of these parameters, to our measurements. 

Again, we use the fitting functions of Smith et al. (2002) to 
produce dark matter power spectra Ps{k,r) for values of as rang- 
ing from 0.1 to 1.5 and fim from 0.1 to 1.0, exploring all values 



on a grid with 0.01 spacing in these parameters. We choose to fix 
ilm + QA ~ 1- For all our theoretical power spectra and correlation 
functions, we use Ho = lOO/i = 68kms"^Mpc"^ for the value 
of the Hubble constant and = 1 for the initial slope of the power 
spectrum of density fluctuations. For each set of parameters, we 
calculate the corresponding shear power spectrum using equation 
(17). Note that we input our composite redshift distribution (Fig. 8) 
into the calculation for the shear power spectra. Having calculated 
the shear power spectra we obtain the correlation functions Ci,2 
using equation (25). 

We fit the above models to our correlation function and shear 
power spectrum data (where the A901 supercluster field has been 
excluded in both cases) using a fitting procedure (c.f. Bacon 
et al. 2002). We order our correlation function measurements as a 
vector d = {Ci{0n), C2{0n)}, where Ci,2(0n) is the mean cor- 
relation function (averaged over all sections) for a given angular 
separation (see Section 5.1). Similarly, we order our theoretical 
correlation functions as vectors x(a) with the same format, where 
a — {as, ftm)- 

We use the log-likelihood estimator 



x' = [d-x(a)]^V-^[d-x(a)], 

where 

V = (dd') 



(52) 



(53) 



is the covariance matrix of our correlation function measurements, 
which we measure directly from our data using equation (33). This 
estimator is valid in the case of Gaussian errors, which we achieve 
since we have averaged over many sections. We calculate for 
our fine grid of {as, Jim) theoretical correlation functions, and find 
the minimum and confidence intervals. 

The fitting for the shear power spectrum measurements is 
done in exactly the same way as for the correlation function analy- 
sis. 



7.1.1 Correlation function results 

For the correlation function analysis, we find a best fit 
as{nm/0.3)°-^^ = 0.75, with reduced of 1.21. The la error 
bar on this value, for a 1 parameter fit, is given by the boundary 
Ax^ = 1-0, which occurs here at as — 0.81 and 0.69. To date, 
one of the biggest sources of error in cosmic shear measurements 
of cosmological parameters has been the uncertainty in the median 
redshift of the source galaxies. The cosmic shear signal scales as 
a?i cc as'^z^'^ (see e.g. BRE). So the uncertainty in the median 
redshift contributes to the error in as as 



Sa's \ _ / Sas\' 



0.64' 



\ as / 



(54) 



Using our estimate of the median redshift of the COMBO-17 galax- 
ies, Zm ~ 0.85 ± 0.05 (see Section 4.2), we have investigated the 
extra uncertainty introduced into our measurement of erg by equa- 
tion (54). We find the additional error introduced due to the me- 
dian redshift uncertainty to be 0.01. Thus we obtain a measurement 
signal-to-noise of ~ 9 for the amplitude of the power spectrum, 
with a measurement of the amplitude 



as{n^/o.sy 



0.75_ 



(55) 
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Figure 19. The likelihood surface of ag and Qm from COMBO-17 as 
calculated using our correlation function measurements (lighter contour) 
and as calculated using our shear power spectrum measurements (dark con- 
tours). For the power spectra measurements we plot the 1 and 2(t contours. 
For the correlation functions, we plot only the Icr contour for clarity. 

7.1.2 Shear power spectrum results 

The results for the shear power spectrum analysis yield a shghtly 
lower amplitude for the power spectrum normalisation. A good fit 
to the parameter constraints as calculated using our shear power 
spectrum measurements is (jg(Slm/0.3)'' '*^ — 0.72, with a re- 
duced of 2.21. The additional error due to the uncertainty in Zs 
is again 0.01. Our final measurement for the normalisation of the 
mass power spectrum, as calculated from our shear power spectrum 
measurements is 

a8(n„/0.3)°-*^ = 0.72i°;°9- (56) 

These constraints, as calculated using both of our methods are 
shown as contours in the erg — flm plane in Fig. 19. Note that in the 
case of the correlation function constraints, values of Qm ^ 0.63 
are excluded at the Icr level. 

7.2 Including the actual redshift distribution 

In order to demonstrate the power of accurate redshift informa- 
tion we have also measured parameters from correlation functions 
which have been calculated using only galaxies with reliable red- 
shift estimates. In this case, for our theoretical curves, we input 
the actual redshift distribution of the lensed source galaxies (i.e. 
the dotted curve in Fig. 8). For the purposes of this demonstra- 
tive calculation, we have included the A901/2 supercluster field. 
We do this simply to increase the number of galaxies - at present, 
photometric redshifts are only available for the CDFS, Sll and 
A901 fields and there were too few galaxies with reliable redshifts 
in CDFS and Sll alone (19143 galaxies) to constrain parameters 
significantly. Adding the 12506 galaxies with accurate redshifts in 
the A901 field to this sample enabled us to obtain reasonable con- 
straints on Qrn and erg. The result, therefore, will clearly be biased 
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Figure 20. The likelihood surface of erg and Qrn, from COMBO-17 using 
correlation functions measured only from galaxies with accurate redshift es- 
timates. This subset of galaxies only includes galaxies in the CDFS, Sll and 
A901 fields with R-band magnitudes < 24.0. The Icr contour is plotted as 
a full line. For comparison, we also plot the Icr confidence region (dashed 
contour) obtained from the same correlation functions excluding the red- 
shift information. For this calculation the only redshift information we as- 
sume is an estimate of the median redshift of the survey (Zm = 0.6 it 0.2). 

by the A901/2 field and cannot be taken as a measure of the power 
spectrum amplitude. However, it does dramatically demonstrate the 
increase in accuracy attainable with accurate photometric informa- 
tion. 

The results are shown in Fig. 20. For comparison, we also plot 
the constraints obtained from the same set of correlation function 
measurements where we have assumed a median redshift and un- 
certainty of Zm = 0.6 ± 0.2. Here, Zm = 0.6 is the measured 
median redshift of galaxies with reliable assigned redshifts in the 
COMBO-17 survey and Az = 0.2 is typical of the estimated un- 
certainty in the median redshift of cosmic shear surveys to date. 
In agreement with the shear power spectrum analysis (see section 
6.3), we see the effect of the A901 supercluster pushing the best-fit 
as value up to ~ 1 for a matter density of ilm ~ 0.3. Once again, 
we emphasise that we are not presenting this as a measurement 
of as - we simply wish to demonstrate the dramatic improvement 
in parameter constraints obtainable with the inclusion of accurate 
photometric information. Clearly for accurate parameter estimates, 
a greater understanding of the source redshift distribution must be 
a priority for cosmic shear surveys in the future. 

7.3 Combination with the 2dFGRS and CMB experiments 

We can combine the confidence region given by equation (56) or by 
Fig. 19, with parameter estimations from other sources such as the 
2dF Galaxy Redshift Survey (2dFGRS, Percival et al. 2002) and 
the various CMB experiments (Lewis & Bridle 2002 & references 
therein). 

In order to combine our measurements with those from the 
2dFGRS and the CMB data, we first impose some priors on the 
various datasets. Firstly, we consider only flat Universe models 



© 2003 RAS, MNRAS 341. 100-1 18 



18 M.L. Brown et al. 



LO 



b 




Figure 21. The likelihood surface of erg and Qm from combining the 
COMBO- 17 dataset with the 2dFGRS and the latest CMB constraints. The 
dark thin solid contours are the constraints obtained from our shear power 
spectrum analysis described in the previous section. The lighter vertical 
contours are the constraints on Q.m obtained from the 2dFGRS where we 
have applied the priors described in the text to the 2dF data. The dashed set 
of contours are the constraints from a compilation of six CMB experiments 
(see text for details) where we have assumed an optical depth to reioniza- 
tion of r = 0.10. The dark heavy contours are the 1,2 and 3(T combined 
constraints from the three methods. The best fit values of Vim and ""8 are 
also indicated. 



(Qa + f2m = 1) as compelling evidence for a flat Universe has al- 
ready been found from the CMB data (e.g. de Bemardis et al. 2002; 
Pryke et al. 2002; Lewis & Bridle 2002). For the slope of the power 
spectrum of the primordial density perturbations, we take Us = 1 
and for the baryon fraction ft = Jit /Jim, we take = 0.16. 
These values are consistent with those measured from joint analy- 
ses of CMB and Large Scale Structure (LSS) data (e.g. Efstathiou et 
al. 2001; Wang, Tegmark & Zaldarriaga 2002; Percival et al. 2002; 
Knox, Christensen & Skordis 2002). To display the additional con- 
straints in the — Jim plane, we also need to choose a value for 
the Hubble parameter, h. We have taken the value h — 0.68, found 
from a joint analysis of the CMB and 2dF data for a flat Universe 
and for a scalar-only spectral index, (Percival et al. 2002; Efs- 
tathiou et al. 2002). Note that this is the same value of h which we 
adopted when making our model predictions in the previous sec- 
tion. Finally, we note that this value is consistent with the HST Key 
Project within the quoted errors (Freedman et al. 2001). 

Having set these various parameters, in the case of the 2dF- 
GRS, for each of our grid points in the (Jim, eg) plane, we sim- 
ply add the corresponding likelihood values from Percival et al. 
(2002) to the likelihoods already obtained in the previous section 
from our shear power spectrum data. For the CMB constraints, we 
have adopted the results of Lewis & Bridle (2002) who provide the 
following joint constraint on Jim and erg as found from a compila- 
tion of six CMB experiments: 

u = Jlmft^-^^(crge~^/0.7)"'' *'' = 0.115 ± 0.0047. (57) 



Here, r is the optical depth to reionization. For each point on our 
(Jim, o"g) grid, we calculate likelihoods for the CMB data as 

XCMB - -^^2 (58) 

with V = 0.115, Siy = 0.0047 and h = 0.68. Finally, we add the 
likelihoods from the various data to get our final measurements of 
Jim and (jg. The results of this combination are shown in Figure 21 
for an optical depth of r = 0.10. Note that the 2dFGRS data we 
have used for this estimation constrains Jim only and so the con- 
straints shown on as come wholly from the cosmic shear and CMB 
measurements. We measure best-fit values of Jim = 0.271q oi and 
as = 0.731q o3 from the combined data. 

We have also investigated the effect a non-zero r has on these 
measurements. We find that increasing/decreasing r has the effect 
of increasing/decreasing the slope of the CMB constraints some- 
what, but the resulting combined constraints are not altered drasti- 
cally. For example, for an optical depth of r = O.OO, the best-fit 
Jim and as values change to 0.29 and 0.72 respectively, while in- 
creasing r to 0.25 changes the best-fit values to Jim ~ 0.25 and 
erg = 0.77. 



8 CONCLUSIONS 

We have presented an analysis of the cosmic shear signal in the 
COMBO- 17 survey, a 1 .25 square degree survey with excellent data 
quality, including photometric redshifts for ~ 38% of our galaxy 
sample for three of the five 0.5° x 0.5° fields. Our measurements 
follow a process of careful data reduction, assessment of the level 
of telescope-induced shear, and correction for PSF anisotropy and 
circularization of galaxies. In this fashion we reduce the residual 
systematic effects to a level < 5% of the cosmic shear signal. 

We have measured shear correlation functions and variance- 
weighted shear cell-variance for our galaxy sample, detecting the 
cosmic shear signal at the 5cr level. On scales ^ 1 arcmin, where 
we are confident of the cosmological origin of the signal, we find 
a somewhat lower amplitude than other cosmic shear surveys have 
measured to date. 

In contrast to previous cosmic shear surveys, we have sub- 
stantially reduced the usual uncertainty introduced into the inter- 
pretation of cosmic shear measurements due to a lack of knowl- 
edge of the redshift distribution. We have done this by estimating 
from the data the median redshift of the lensed source galaxies to 
be z,n = 0.85 ± 0.05. 

We have used our catalogues of galaxy positions and shear 
estimates to apply a maximum likelihood analysis, obtaining the 
shear power spectrum for each COMBO-17 field and for the entire 
survey. We have measured the K-field power spectrum expected 
to arise from gravitational shear, from £ = 400 to 10*, and have 
shown that the systematically-induced /3-field power spectrum is 
substantially below our shear signal throughout this £ range. We 
find that the power spectrum for a supercluster field has a signif- 
icantly higher amplitude than that for random fields, as expected, 
and we have therefore excluded this field from our final optimally 
combined result. Simulations of the maximum likelihood proce- 
dure demonstrate that our method is unbiased to within ~ 10%. 
We have measured the covariance matrix of our shear power spec- 
trum band powers directly from the data and find our band powers 
to be essentially independent measured of power. 

We have also investigated the use of correlation function es- 
timators (e.g. SvWKM) for the shear power spectrum and find, in 
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agreement with SvWKM, that their success is dependent on the cor- 
relation functions being measured over a large range of scales. Such 
estimators are therefore difficult to apply to datasets composed of 
small fields. 

We have used our shear correlation function and power spec- 
trum measurements to estimate constraints on the cosmological pa- 
rameters f2m and as. We have made fits of our measurements 
to predicted functions for various values of these cosmological pa- 
rameters, and we measure the amplitude of the mass power spec- 
trum, from our shear power spectrum results to be crs (^j^) = 
0.72 ± 0.09 with 0.10 <Q.m< 0.63, including uncertainties due 
to statistical noise, sample variance, and redshift distribution. 

We have demonstrated the potential of including photometric 
redshift information by constraining these cosmological parameters 
with only those galaxies with assigned redshifts. For this calcula- 
tion, we have input the actual measured redshift distribution into 
the theoretical predictions and have compared the resulting con- 
straints with those that would be obtained with only an estimate 
of the median redshift of the survey. We thus demonstrate the dra- 
matic increase in precision obtained by including accurate redshift 
information for the source galaxy distribution. 

We have combined our constraints with those from the 2dF 
Galaxy Redshift Survey, and the latest CMB experiments, finding 
a power spectrum normalisation of as = 0.731q q3 and a matter 
densityof n„ = 0.27lo:^i. 

These results (for both the lensing alone, and for the com- 
bined constraints) are lower than previous constraints found from 
cosmic shear surveys (e.g. Bacon et al. 2002, Refregier et al. 2002, 
Van Waerbeke et al. 2002, Hoekstra et al. 2002) which have, until 
now found power spectnmi normalisations of as ~ 0.85 — 1.0. 
However, our measurements are much more in agreement with re- 
cent cluster abundance estimates of the power spectrum normal- 
isation (Borgani et al 2001; Seljak 2002; Reiprich & Bohringer 
2002; Viana, Nichol & Liddle 2002, Allen et al. 2002; Schuecker 
et al. 2003, Pierpaoli et al. 2002) and those found from combin- 
ing constraints from the 2dFGRS and CMB results (e.g. Lahav et 
al. 2001, Melchiorri & Silk 2002). Finally, we note that our results 
agree within the 2a level with recently aimounced results from the 
WMAP satellite (Spergel et al. 2003). 

Since this paper was submitted, two other groups have pre- 
sented cosmic shear results yielding somewhat lower values of 
an - Jarvis et al. 2003 measure the shear signal from the 75 sq. 
degree CTIO survey and find a power spectrum normalisation of 
as = 0.71 while Hamana et al. 2002 have used the Suprime-Cam 
instrument on the Subaru telescope to measure the cosmic shear 
signal from 2.1 sq. degs. of deep R band data, with which they 
constrain the power spectrum normalisation to be as = 0.69. The 
results from both these studies are in excellent agreement with the 
results presented here. 

We have tested our dataset extensively for systematic effects 
by measuring galaxy cross-correlations, star-galaxy correlations 
and by decomposing the cosmic shear signal into its constituent 
curl and curl-free modes. However, we have foimd our dataset to 
be largely free of any systematic effect that could account for the 
discrepancy found between our parameter constraints and those ob- 
tained from previous weak lensing studies. One factor, which must 
be important, is the increased understanding of the source redshift 
distribution which we have gained from the 17-band photometric 
information contained in the COMBO- 17 survey. With the increase 
in area of cosmic shear surveys in the future, this type of redshift 
information is likely to become vital as the uncertainties due to the 



redshift distribution become dominant over other sources of error 
such as shot noise and sampling variance. 
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